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ABSTRACT 

We present a detailed study of the nuclear structure of the highly polarized broad absorption line 
quasar, Mrk 231, through a multiwavelength campaign of Chandra observations, optical spectroscopy, 
optical spectropolarimctry, and imaging polarimetry. This campaign was designed to extend the 40 ks 
Chandra study of Gallagher et al. and the optical and UV spectropolarimetric study of Smith et al. 
to probe variability on multiple timcscalcs. As direct emission from the nucleus is heavily obscured 
at optical through X-ray wavelengths, the detailed study of scattered emission has lead to insights 
into the stratified and complex central region of this active galaxy. Though significant continuum 
variability is not detected in any of the three new 40 ks Chandra observations, we investigate Fe Ka 
emission features in the individual and combined spectra. Comparing echelle spectra of the Na ID 
absorption lines with the literature, we show that one system disappeared in 2000 only to reappear in 
later epochs. Notably, we detect a large decrease in polarization across the entire optical spectrum of 
Mrk 231 from 1995 to 2003. Though the polarization fraction fell, e.g., from 4% to 3% at A obs ~ 6000A, 
the polarization position angle spectrum remained unchanged. The optical polarization behavior is 
consistent with a decrease in the flux scattered by circumnuclear material on spatial scales where the 
broad-line region is resolved by the scattering material. Ultraviolet imaging polarimetry of Mrk 231 
by the Hubble Space Telescope sets an upper limit on the distance between the active nucleus and the 
scattering regions of <20pc. 

Subject headings: galaxies: individual (Markarian 231) — galaxies: nuclei — galaxies: Seyfert - 
polarization — ultraviolet: galaxies — X-rays: galaxies — quasars: absorption 
lines 



1. INTRODUCTION 

The variability of emission from active galactic nu- 
clei (AGNs) has long been employed in efforts to dis- 
cern the structure and physics of material near the cen- 
tral energy source. The dramatic changes in brightness 
and polarization of the continua of radio-loud objects, 
particularly at short wavelengths, played an important 
role in developing the popular model of a relativistically 
beamed jet emitting by t he synchrotron mechanism (e.g., 
IBlandford fc^eelfl97l . The effects of more isotropic 
changes in luminosity have been exploited by the "re- 
verberation mapping" technique to infer the size, thick- 
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ness, and ionization properties of the broad emission-line 
region, and more recently to identify systematic trends 
in nu clear parameters between various classes of AGNs 
fe.g..lNetzer fc PetersonllT997h . 

Evidence for changes of the structure of the nuclei of 
active galaxies can be found in the appearance and dis- 
appearance of shoulders and other features in the pro- 
files of broad emission lines, modifications of the B aimer 
decrements, and in certain cases in the replacement 
of a normally centrall y peaked profile b y one that is 
strongly double-peaked l|Ulrich et alll997l and references 
therein). Some of these effects appear to be induced by 
luminosity changes of the central continuum source, but 
others require mass motions or the redirection of an il- 
luminating beam, also known as the 'searchlight' effect. 
At larger radii, variations in the properties of absorbing 
gas as indicated by changes in UV absorption lines have 
been detected on timcscales of ~1 yr in several Seyfert 
galaxies, with explanations including both luminosity- 
induced ionization changes and the motion of absorbers 
across our line of sight fe.g.. lGabel et alll2003l and refer- 
ences therein). Similarly, large column density variations 
have been observed in X-ray studies of absorbed AGNs 
on timescales of mont hs to years fe.g.. lRisaliti et alJ2002t 
Gallaghe rTTaT]l200l . 

The luminous Type 1 (broad emission line) AGN 
Markarian 231 is well known to exhibit no less than 
three displaced o ptical absorption-line sy stems in Na I, 
He I, and Ca II l|Boksenberg et al.lll977l(l and a fourth 
appe ared in observations made in 1988 ( Bor oson et alJ 
I1991|) . Further variability in the strengths of the ab- 
sorption lines has been reported in subsequent years 



2 



Gallagher et al. 



QKollatschnv et all Il992t iForster et aD 119951 hereafter 
FRM) and Mg IIA2800 ab sorption associated with one 
of the systems reported bv lSmith et al.1 l)1995t hereafter 
SSAA). This observation, coupled with the source's high 
intrinsic luminosity, qualified Mrk 231 as a low- ionization 
broad absorption-line (LoBAL) quas ar. LoBAL quasars 
are notoriously weak X -ray emitters l)Green et al.1 1200 ll 
iGallaeher et al.1 1200 2a ): therefore the proximity (z = 
0.0422) of Mrk 231 makes it by far the most easily stud- 
ied example of this class in X-rays. 

The fact that most of the ~4 x 10 12 luminos- 
ity o f Mrk 231 appears in the infrared l|Rieke fc Low! 
1972) was early evidence that the central energy source is 
heavily obscured along most sightlines. This is strongly 
supported by the paucity of emission from the narrow- 
line region (specificall y, the absence of [O III] flux; 
Hamilton fc KeeH IT987). and the shape and luminosity 
of the 2-10 keV X-ray spectrum imply that the nu- 
cleus is only seen via reflection and/or scattering in that 



band i|Malonev fc Revnoldsl2 000: Gallagh er et al.l2002bl 
hereafter G02). An analysis of 15-60 keV Beppo-SAX 
Phosw ich Detector System (PDS) data bv iBraito et al.l 
(2004) confirmed this interpretation with the detection, 
for the first time, of the direct X-ray continuum absorbed 
by an intrinsic column density of ATh~ 2 x 10 24 cm~ 2 . 
Absorption of this magnitude completely blocks the di- 
rect 2-10 keV continuum, which has an absorption- 
corre cted luminosity L 2 -io ~ 10 44 ergs s _1 l)Braito et al.l 



2004). The observed, nonthermal 2-8 keV continuum, 
with Lx ~ 2 x 10 42 ergs s _1 , must therefore be repro- 
cessed emission. However, the 2-8 keV flux variability 
on timescales of hours detected by G02 with Chandra 
indicates that the X-ray scattering medium is compact, 
with spatial scale of ~ 10 15 cm. Given the proximity 
of this gas to the central continuum source, this scat- 
tering medium is likely to be highly ionized. Obscu- 
ration and the attendant scattering of nuclear light by 
pervasive dust clouds at much larger radii, R > 1 pc, 
give rise to polarization of the optical- UV light, and 
polarization fractions that reach 15% in the near-UV 
have made the object a popular target for polarimet- 
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The sensitivity of both the X-ray spectrum and the 
optical polarization of Mrk 231 to the geometry of scat- 
tering on different physical scales, and the insight that 
variations in these diagnostics might lend to our knowl- 
edge of the structure of the nuclear regions of a lumi- 
nous Type 1 AGN, have led us to carry out a program 
of high-quality, ground-based optical spectropolarimetry 
and 0.5—8 k eV spectroscopy with the Chandra X-ray 
Observatory l)Weisskopf et alJ 12002) . In this paper we 
present the results of that study as well as optical echelle 
spectroscopy and UV polarimetry from different epochs 
that contribute additional information on the properties 
of the nucleus. 

The cosmology assumed throughout has Hq = 
70 km s" 1 Mpc" 1 ', Q M = 0.3, and fi A = 0.7. For Mrk 231 
at z — 0.0422, 1 arcsecond corresponds to a physical scale 
of 833 pc. 



2. X-RAY OBSERVATIONS AND ANALYSIS 

Mrk 231 has been observed 4 time s with Chandra, each 
time at the aimpoint of the ACIS ijGarmire 
back-illuminated S3 detector. The first observation of 
2000 Oct 19 (observation 1031) was in full-frame mode. 
The subsequent three observations of 2003 Feb (obser- 
vations 4028-30) were performed in 1/2 subarray mode 
in order to reduce the ~ 7% pile-up fraction of the first 
observation (G02). Each ACIS observation was between 
36-39 ks long; Table ^ presents the log of the observa- 
tions, including dates, exposure times, and soft and hard- 
band count rates. The At of ~lweek between observa- 
tions 4028, 4029, and 4030 was chosen to probe times 
with a scale intermediate between the hours and years 
time scales of previous observations. Physically, this cor- 
responds to structure on scales of ~ 10 16 c m, the scale of 
the th ick, highly ionized gas postulated bv lMurrav et alJ 
( 1995) to shield a radiatively driven BAL wind from be- 
coming overionized by the ionizing X-ray continuum. 

The Chandra data were processed with the standard 
Chandra X-ray Center (CXC) pipeline from which the 
level 1 events file was used. The data were filtered on 
good time intervals. The standard pipeline processing in- 
troduces a 0'.'5 software randomization into the event po- 
sitions, which degrad es the spatial resolution by ~ 12% 
ijChartas et al.l 12002^1 . This randomization is done to 
avoid aliasing effects noticeable in observations with ex- 
posure times of less than 2 ks. Since spatial resolution 
was important for our analysis and the ACIS observa- 
tions were significantly longer than 2 ks, we removed 
this randomization from the ACIS observation by recal- 
culating the event positions in sky coordinates using the 
CIAO tool, acis-process_events. The ACIS-S3 imaging 
data were corrected for radiation damage sustained by 
the CCDs during the early mont hs of Chandra op e ration s 
using the software described in iTownslev et all (2002). 
This procedure mitigates the effects of charge transfer 
inefficiency (CTI) that introduces positional dependence 
of the gain, quantum efficiency, and grade distribution. 
CTI-corrected ACIS-S3 data allow a single redistribution 
matrix file (rmf) to be used across the whole CCD. Af- 
ter CTI-correction, the data were filtered to keep "good" 
time intervals and ASCA grades (grades 0, 2, 3, 4, and 6). 
Filtering on "status" can remove X-ray events from the 
nucleus because X-rays can be inappropriately tagged as 
cosmic-ray afterglows. Therefore, the values in the sta- 
tus column for cosmic-ray afterglows were ignored for the 
final filter on "good" status. 

2.1. Variability Analysis 

Given the significant hard-band X-ray variability de- 
tected in the first Chandra observation, we first searched 
for count-rate changes both within and between the ob- 
servations. The time and energy of each event within a 
3"-radius source circle centered on the peak of the X-ray 
emission were extracted from the reduced events file for 
each observation. The unbinned full (0.5-8.0 keV), soft 
(0.5-2.0 keV), medium (2.0-5.0 keV), hard (2.0-8.0 keV), 
and very hard (5.0-8.0 keV) band event arrival times for 
each observation were evaluated with the non-parametric 
Kolmogorov-Smirnov test to search for variability within 
each observation. For the new observations, no signifi- 
cant variability was detected in any band. The original 
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data from 2000 (observation 1031) were downloaded from 
the public archive, rereduced with CIAO 3.0 as outlined 
above, and retested. As presented in G02, we recovered 
significant variability in the medium and hard bands at 
the > 99% confidence level; thi s result was also co nfirmed 
by the independent analysis of iPtak et alJ l|2003f) . 

To check for variability between observations, the 
background-subtracted nuclear count rates were com- 
pared in the five bands. The chosen background region 
was a source-free annulus with an inner radius of 40" and 
an outer radius of 50", well outside of the galaxy, cen- 
tered on the nucleus. Figure Q] presents these data. The 
full-band count rate is remarkably stable, but there ap- 
pears to be a slight shift over the course of observations, 
with more hard-band counts relative to soft-band counts 
in the last observation than in the preceeding ones. This 
trend is slight however, and we investigate any possible 
spectral hardening further with spectral analysis. 

2.2. Nuclear Spectroscopy 

The nuclear spectra were extracted from each reduced 
Chandra dataset independently. The source and back- 
ground regions were the same as those used in the vari- 
ability analysis. The auxiliary response files (arfs) and 
rmfs were generated with the CIAO 3.0 tools mkarf and 
mkrmf, respectively. 

The final, preferred model from G02 for the fit to the 
Chandra spectrum included a thermal component from 
the nuclear starburst as well as three absorbed power-law 
components with different normalizations and intrinsic 
column densities (see Figure 9d of G02). Given the ob- 
served 2-8 keV luminosity of the source (~ 100 times 
fainter than expected based on both the dereddened UV 
luminosity and the > 10 keV X-ray spectrum) and the 
known complexity of the scattering geometry derived 
from optical polarization studies (SSAA), these power- 
law continua were interpreted as direct continuum X-rays 
scattered into the line of sight by distinct structures, and 
therefore passing through different amounts of absorp- 
tion. The scattering was inferred to be by electrons in 
highly ionized plasma, thus preserving the shape of the 
intrinsic spectrum. In G02, the photon index of each 
power-law model was fixed to T — 2.1, based on the cor- 
relation of H/3 FWHM and T found bv lReeves fc Turned 
(2000) from ASCA spectral fitting of quasars. To evalu- 
ate this model in the context of the new Chandra data, 
we examined the combined, average spectrum from the 
three new observations to increase the signal-to-noise ra- 
tio. With these data, the power-law photon index can be 
fit independently. We do not incorporate the data from 
the first Chandra epoch because they were mildly piled- 
up and taken in a different mode (full array rather than 
1/2 subarray). 

To create an average nuclear spectrum, the three un- 
binned spectra and background spectra were each added 
together using the FTOOL mathpha. An average arf 
and rmf were created by weighting by the 0.5-8.0 keV 
counts using addarf and addrmf, respectively. The com- 
bined spectrum (with a total exposure of 114.3 ks and 
4300 counts) w as then fit usin g the X-ray spectral-fitting 
tool, XSPEC (|Arnaudlll99ej) . Three absorbed power- 
law models plus a Raymond-Smith plasma were ap- 
plied to the combined spectrum; T was tied for each 
power law. An acceptable fit {y^/v = 174.8/153) was 



found for this model with T = 1.923o7o- However, 
given the relatively featureless spectrum, this solution 
is not unique, and much lower values of T also provide 
statistically acceptable fits. As T ~ 1.9 is typical of 
radio-quiet quasars and more plausible in the context 
of this model, we fixed T = 1.92 to determine the er- 
rors in the remaining parameters. The column densities 
(assumed to be neutral) for the three power laws are 
iVH=(0.19±g;g|, 4.6±g;|§, 45±^) X 10 22 cm" 2 , similar to 
the values listed in Table 2 of G02. This model does 
not include a component corresponding to the highly 
absorbed (Vh~ 2 x 10 24 cm~ 2 ) d irect X-ray contin- 
uum studied by IBraito et alJ l|2004l) because this com- 
ponent does not contribute significant counts to the 
Chandra spectrum. The best-fitting temperature for the 
Raymond-Smith plasma is kT = 0.89±g;$j keV. Though 
the spectral bump near 1 keV (likely from spectrally 
unresolved Fe L emission lines) is fit by this model, it 
is also possible that t here is a contribution from pho- 
toionized plasma, and iTurner fe Kraemerl 1)20031 here- 
after TK) suggested that an emission feature seen at 
~ 0.45 keV in the XMM-Newton spectrum might be a 
C IV radiative recombination line. 

TK detected a broad, FeKa emission line with an 
equivalent width EW~ 450 eV, a rest- frame energy 
E = 6.49 ± 0.12 keV, and a = 0.17 ± 0.12 keV in the 
17.3 ks XMM-Newton EPIC-pn spectrum of Mrk 231. 
Such a strong, broad feature is not evident in this 0.5- 
8.0 keV combined spectrum, and adding a broad line does 
not make a significant improvement in the spectral fit 
(A% 2 = —2.8 for three fewer degrees of freedom). Adding 
a single, unresolved line in the Fe region has a best-fitting 
energy of 6.641°'^ keV and gives A\ 2 — —5.5 for two 
fewer degrees of freedom. 

To focus further on the nuclear spectrum in general 
and the Fe Ka region in particular, we follow the exam- 
ple of TK by ignoring the data below 3.0 keV where star- 
burst emission is significant and fitting the 3.0-8.0 keV 
counts with a single power-law model for the continuum. 
From G02, this is also the region of the spectrum that 
showed significant variability in observation 1031 of 2000 
Oct. To account for the FeKa edge absorption expected 
from multiple absorbers (not otherwise included in this 
model), we also add an edge with the rest-frame energy 
fixed at 7.12 keV, the ionization energy of neutral Fe. 
This two component model is sufficient to fit the data, 
with x 2 /^=69.0/69, and adding a broad line does not 
significantly improve the fit (A% 2 = —3.2 for three fewer 
degrees of freedom). Convolving the best-fitting XMM- 
Newton model with the Chandra response and allowing 
F, the continuum normalization, and r max (the maxi- 
mum optical depth of the ionization edge) to vary while 
keeping the Fe Ka line parameters fixed does not provide 
a better fit (\ 2 jv = 72.8/69) than the power-law plus 
edge model alone. Freeing the continuum parameters 
is necessary to account for both cross-calibration uncer- 
tainties and potential source variability. A strong nega- 
tive residual at rest-frame ~6.45 keV (approximately the 
peak energy of the XMM-Newton-detected emission line) 
suggests that Chandra is resolving the broad line seen by 
TK into two narrow lines. At 6 keV, the XMM-Newton 
EPIC-pn has an effective area about 4 times larger than 
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that of ACIS-S3. 9 Our 114 ks observation is about 6.6 
times longer than the XMM-Newton observation, and 
therefore the Chandra data have better photon statis- 
tics. Both the Chandra and XMM-Newton data can 
be successfully fit with two unresolved (a = 0.01 keV) 
gaussian lines, with consistent line energies of ~6.3 and 
~6.6 keV. In the XMM-Newton data, with fewer counts 
in the observed-frame 5.8-6.7 keV band (~ 140 versus 
~ 210 for the combined Chandra spectrum), two gaus- 
sians are not preferred over a single, broad gaussian. Fig- 
ure shows the combined data and preferred spectral 
model, and the parameters from the preferred spectral 
models for both the Chandra and XMM-Newton data 
are presented in Tabled 

From visual inspection of a merged 5.8-6.7 keV (rest- 
frame 6-7 keV) image, a 3"-radius source cell encom- 
passes all of the counts in this band. That is, there is no 
evidence from the Chandra data of significant Fe-band 
emission coming from the host galaxy. Therefore, any 
discrepancies between the observation presented in TK 
and the Chandra data do not arise from the different 
beam sizes of Chandra and XMM-Newton. 

We investigate further the possibility of nuclear spec- 
tral variability between the 5 observations (four Chandra 
and one XMM-Newton) . Again, each spectrum was ini- 
tially fit jointly with a power law model and an absorp- 
tion edge fixed to 7.12 keV. Both T and r max were tied 
for the fit, and the power-law normalizations were free 
to vary for each observation. To look for iron emission, 
the data in the iron band (rest-frame 6-7 keV, observed- 
frame 5.8-6.7 keV) were ignored for this continuum fit. 
The Chandra residuals from this joint-spectral fitting are 
shown in Figure |21 For the final observation, there is no 
strong evidence of iron emission. For the first observa- 
tion, a series of positive residuals in the iron band is 
suggestive of a broad emission feature, as found by TK. 
In addition, for the middle two observations, Feb 03 and 
Feb 11, there are consecutive positive residuals at the 
1—2(7 level in the iron-band region. To examine these 
features more closely, we fit the spectra individually. 

For each observation, the power-law model with an 
edge at 7.12 keV was applied to the data. The values of T, 
T mH , and the power-law normalization were free to vary. 
In each case, the fit was acceptable with x 2 /^=29.8/30, 
29.1/31, 30.3/33, and 13.1/31 for observations 1031, and 
4028-30, respectively. Next, a narrow iron line was added 
to each model as a redshifted gaussian with the energy 
free to vary between rest-frame 6-7 keV. The value for 
a was fixed to 0.01 keV, unresolved at ACIS CCD res- 
olution. In chronological order, the new values for % 2 
were 26.8, 22.4, 21.7, and 13.1 for two fewer degrees of 
freedom. Of these, only the middle two (observations 
4028 and 4029) indicate a significant improvement in the 
model fit with A% 2 = —6.7 and -8.6, respectively. Ac- 
cording to the F-test, the addition of the gaussian com- 
ponent improved the individual models at the > 97% 
and > 99% confidence level for the data from observa- 
tion 4028 and 4029, respectively. For observation 1031, a 
broad line with a free to vary was also tested. For three 
fewer degrees of freedom from the power-law plus edge 
model, Ax 2 = —6.6, indicating an improvement at the 
~ 92% level over the continuum model alone. Formally, 

9 http:/ /heasarc. gsfc.nasa.gov/docs/xmm/about_why.html 



these confidence levels should be lower, as they apply to 
a single observation (rather than the four we are present- 
ing), but they are still useful as a relative gauge of the 
improvement. The values and errors for the parameters 
of the preferred model for each spectrum are listed in 
Table □ 

In Figure the contour plots for the line flux versus 
line energy for 4028 and 4029, the observations with two 
possible narrow lines, are plotted. As is evident from the 
figure, the 68% and 90% contours do not overlap; the 
energies are not consistent. The fluxes and EWs found 
for both observations are however consistent. 

2.3. X-ray Properties of the Galaxy 

A merged event file of the three new observations was 
created using the CIAO script merge-all. This script re- 
projects the positions of events from all images into a 
common frame so that they can be combined. The im- 
ages were then adaptively smoothed using the CIAO tool 
csmooth. The X-ray emission from the galaxy based on 
the 40 ks 1031 observation was described in detail by 
G02, and so we only comment briefly on the new infor- 
mation obtained from these additional data. 

From the full-band image (Figure EJ, the deeper data 
allow us to see ring structure around the nucleus. A 
radial surface-brightness profile of the soft-band (0.35- 
2.0 keV) counts centered on the nucleus (Figure[|JJ| reveals 
an enhancement of the X-ray emission at R rj 2'.' 6. This 
radius is approximately halfway between th e two b right 
optic al 'shells', SI and S2 (as labeled by iLfpari et alJ 
2005), seen in HST images. The overall galactic emis- 
sion is structured and asymmetric, with X-ray enhance- 
ments to the northeast and a deficit to the west of 
the nucleus; both were previously seen in the 1031 im- 
age. The northeast structure is likely from photoion- 
ized (rather than collisionally ionized) gas as it coin- 
cides with an optically bright regi on attributed to [O III] 
emission (|Hutchings fc Nef_ll987h . In the full-band im- 
age, the off-nuclear point-source detected by G02 is no 
longer present, however, a second one (only evident at 
energies < 3 keV) has appeared to the southwest of the 
nucleus. The inferred variability of both sources is con- 
sisten t with an identificatio n as luminous X-ray binaries 
(e.g.. iFabbian o et alJ 12003]) . preferentially found in ac- 
tively star-forming galaxies. 

The three-color image, made by combining 0.35-1.2 
(red), 1.2-3.0 (green), and 3.0-8.0 keV (blue) csmoothed 
images illustrates the energy dependence of these fea- 
tures (Figure |3J. The asymmetry of the hard emission 
surrounding the nucleus causes the image to appear bluer 
to the west of the nucleus. Except for a compact region 
within ~ 4" of the nucleus, the galaxy is not detected at 
energies >3 keV. 

3. OPTICAL OBSERVATIONS AND ANALYSIS 

3.1. Echelle Spectroscopy 

Mrk 231 was observed with the 9. 2 m Hobby-Eberly 
Telescope (HET; IRamsev et aIIIT99l at McDonald Ob- 
servatory using th e Upgraded Fiber O ptic Echelle Spec- 
trograph CUFOE; lHarlow et al.lll996|) in 2000 January 
and April. These observations were made to search for 
relatively short-time scale (weeks to months) variations 
in the prominent Na I D absorption systems that are 
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blueshifted 5000-6000 km s _1 from the systemic velocity 
of the host galaxy. All of the data were acquired using 
a 400 /im (2") fiber with sky subtraction accomplished 
with a second 2" fiber offset from the galaxy. The spectra 
have a resolution of ~0.8 A in the echelle order contain- 
ing the Na I D absorption features. As summarized in 
Tabled the 2000 January observations consisted of three 
200 s exposures of the nucleus on January 6 followed by 
four 900 s exposures on January 13. No significant dif- 
ferences were seen between the two sets of observations 
and the results were averaged. Similarly, the 2000 April 
3 dataset is an average of three 900 s exposures. 

The January and April data were then inspected 
for differences in the positions, widths, and equivalent 
widths (EWs) of the various absorption systems of Na I 
D. The lower spectrum shown in Figured is the differ- 
ence spectrum (Jan — Apr) between the two epochs, after 
normalizing the spectra to the surrounding continuum. 
Again, no significant differences are seen between the two 
epochs, and the average of the two spectra is displayed as 
the upper trace in Figure [7] with individual line systems 
labeled according to the notation of FRM and quantified 
in Table 01. At the resolution and signal-to- noise ratio of 
the HET observations, the doublet splitting is cleanly re- 
solved and subsystems can be recognized in both systems 
I and II. The Di component of system Ii is a tentative 
identification within the main trough of system I at about 
the correct position relative to D 2 . Expected positions o f 
the lines for the system III doublet ijBoroson et al.in~9Ql^l 
are marked, but it is clear that these features were absent 
in early 2000. 

3.2. Ground-Based Imaging and Spectropolarimetry 

Optical spectropolarimetry of Mrk 231 was obtained 
on several nights during 2002—2004 w ith the spectropo- 
larimeter SPOL ijSchmidt et al.lll992j) . This is the same 
instrument that was used for the observations published 
by SSAA from 1992—1993, and all acquisition, calibra- 
tion, and reduction procedures are common to the two 
epochs. As with other dual-beam designs, the instru- 
ment achieves high throughput while accurately can- 
celling sky background and avoiding systematic effects 
to better than 1 part in 1000. Calibration of the polariz- 
ing efficiency of the optics is performed by observations 
of a continuum source through a fully polarizing prism 
each observing run. However, the use of different in- 
terstellar polarization standard stars for the various ob- 
serving runs can result in a dispersion of up to 1° in 
the overall registration of the position angle. The de- 
tector was also upgraded twice in the interval between 
the epochs, providing improvements in sensitivity, focus, 
spectral coverage, and flat-field calibration. The signifi- 
cant changes for polarimetry are an enhancement of the 
signal-to-noise ratio for a given exposure time and an 
improvement in spectral resolution. As summarized in 
the observation log in Table |31 virtually all data from 
the two epochs were gathered at the 2.3 m Bok telescope 
of Steward Observatory using a low-resolution grating 
(AA ~ 17A) with wide spectral coverage, though one 
dataset in each epoch utilized a 1200 line mm" 1 grating 
with a resolution of 8 A to record better the polarization 
profiles of the Ha emission and Na I D absorption fea- 
tures. In addition, observations were made at the MMT 
in 2004 April to investigate the diluting effects of the 



host galaxy, which spans > 1' ^Hamilton fc Kee1lll987j) . 
Unfortunately, the utility of these data was compromised 
by poor (~2" FWHM) seeing at the time. 

We have intercompared the polarimetric results among 
the individual runs of 2002—2004 in both detail and for 
broad wavelength intervals such as those listed in Tabled 
While the effects of variable host galaxy dilution cause 
polarization changes at the level of AP = 0.1 — 0.2% 
(cf. the results for the two widely different aperture sizes 
from consecutive nights in 2003 Feb), we find no sig- 
nificant variability that can be attributed to the AGN 
during the 16 month period. A reinspection of the data 
from 1992—1993 leads to the same conclusion for that 
period. However, there is compelling evidence for a de- 
cline of ~ 30% (AP = 1 — 2% across the bandpass) in the 
degree of polarization during the decade that separates 
the two epochs. At the same time, the position angle 
remained constant. Means for the AA6100 — 6600 contin- 
uum interval are P = 3.41% ±0.04% at = 96?9±0?4 for 
the 1990's versus P = 2.52% ± 0.10% at 9 = 97?8 ± 0?7 
for the 2000's, where the uncertainties are estimated from 
the scatter (rms) of the individual measurements for each 
epoch. The difference is easily seen in the mean polariza- 
tion spectra for the two periods compared in Figure 03 
and in the binned continuum values listed in Table 
The spectropolarimetric observations were obtained un- 
der a variety of conditions and the emphasis was not on 
accurate absolute flux calibration, but based on the gen- 
eral consistency in flux level between the two epochs, any 
systematic change in total brightness can be limited to 
<15%. 

Despite the change in fractional polarization level, the 
essential polarimetric characteristics of Mrk 231 persist 
between the decades. These include a slow rotation in 
position angle 9 over the spectrum, significant changes 
in 9 between the Fe II and Balmer lines versus the con- 
tinuum, a polarized flux spectrum P x F\ that is bluer 
than the total flux spectrum (P\ oc A~ 2 ), and emission- 
line features that are much narrower in polarized versus 
total flux. 

In the interval between the two spectropolarimetric 
epochs, broad-band imaging polarimetry of Mrk 231 
was obtained at optical and near-UV f H3.3|l wavelengths. 
Ground-based B and R observations were made at the 
Bok 2.3 m telescope in late 1996 with SPOL in its imag- 
ing mode (a plane mirror is inserted behind the Wollas- 
ton prism instead of the diffraction grating), and the re- 
sults are summarized in Table |3| The scale of the SPOL 
images is 0'.'52 pixel -1 . The entries for P and 9 in Table|3] 
are for a 3" x 10" aperture centered on the nucleus in the 
51"x51" images. This rectangular aperture, oriented at 
a position angle of 90°, was chosen to mimic the slit and 
east-west extraction apertures used for spectropolarime- 
try. The optical broad-band results for the polarization 
position angle are consistent with both epochs of spec- 
tropolarimetry, and also show that Or — 0b = 2-3°. The 
level of polarization in late 1996 is intermediate between 
the higher P measured in the early 1990s and the lower 
levels observed after 2001. We have searched the litera- 
ture without finding accurate polarimetry of Mrk 231 in 
this interval. Thus, while it is likely that the polariza- 
tion change in Mrk 231 was occuring during this period, 
we do not have data to determine how rapidly P varied, 
though the intermediate level of the data from 1996 sug- 
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gests that the change did not occur on very short time 
scales. 

The effect of starlight from the host galaxy on the po- 
larization measurements of Mrk 231 can be estimated 
using the imaging polarimetry. Figure 1 101 shows P as 
a function of aperture radius for both filter bandpasses. 
The measured polarizations at B and R for apertures 
ranging from — 1"- 10" in diameter readily show the rel- 
ative increase in starlight with aperture size. Extrap- 
olation to r = 0" of third-order polynomial fits to the 
measured polarization indicates that the intrinsic polar- 
ization of the nucleus is 6.45% and 3.96% in the B and R 
bands, respectively. The ratios of star-to- AGN light us- 
ing a 10" aperture are 0.21 (B) and 0.64 (R). Within the 
typical spectroscopic extraction aperture used for spec- 
tropolarimctry (3"xl0"), the ratios are —0.13 (B) and 
-0.52 (R). 

3.3. HST Imaging Polarimetry 

On 1998 November 28 (UT), imaging polarimetry of 
the inner — l"8xl"8 of Mrk 231 was obtained using the 
Faint Object Camera (FOC) aboard the Hubble Space 
Telescope (HST). These observations were made in an 
attempt to resolve and identify the scattering regions 
around the nucleus. The f/96 optics of the FOC were 
used and the detector was read out in its smallest for- 
mat (128x128 pixel) to mitigate saturation problems 
from the bright nucleus. To minimize the contribution 
of host-galaxy starlight, the observation was made in the 
near ultraviolet with the F346M filter (A = 3480 A; 
AAfwhm = 434 A) . This choice of filter also took advan- 
tage of the high polarization of the nucleus (P — 15% at 
A - 3000 A) measured by SSAA. 

Four exposures ranging from 1166 s to 1363 s were 
made with each of three Polaroid filters in the beam 
combined with F346M. About 1.4 hr was spent using 
each of the POL60 and POL120 Polaroid filters, while 
— 1.5 hr total exposure time was used for POLO. In gen- 
eral, care must be taken with FOC polarimetry because 
the throughputs of the polarizer filters are not identical; 
however, in the F346M bandpass they are quite similar. 
In addition, the instrumental polarization introduced by 
off-axis optics inserted during the 1993 Space Shuttle 
HST servicing mission to correct th e telescope's spherica l 
aberration is not well determined ijCapetti et al.l ll996). 
Given the high polarization of Mrk 231, we have ignored 
the 1-3% instrumental polarization in the analysis of 
the FOC data since any resolved polarized regions will 
have even higher polarization than the —15% observed 
in the near UV for the integrated AGN light (SSAA). 
T he polarimetry was r educed in the manner described 
bv lCapetti et all <|1995t) . 

The imaging data do not show any evidence for ex- 
tended structure. Identification of faint extended emis- 
sion very close to the point source is made especially dif- 
ficult because each Polaroid filter yields a distinct point 
spread function (PSF). No PSF calibration images ex- 
ist for the individual filter combinations, and so accu- 
rate subtraction of the PSFs is not possible. However, 
aperture polarimetry of the images supports the conclu- 
sion that at the resolution of the FOC, Mrk 231 is unre- 
solved in the near UV. Circular apertures with radii of 
3 pix (0"043) to 33 pix (0"469) all yield polarizations of 



P ks 15.5% at 9 w 99°. The measurement of P is close 
to the value found at —3000 A in a much larger aperture 
by SSAA, though the measured position angles are dis- 
crepant by nearly 10°. The difference could be caused 
by the instrumental polarization of the FOC. Also, the 
same polarization values are found for the faint emission 
surrounding the core of the PSF. Aperture polarimetry 
using an annulus with an inner radius of 6 pixels from 
the centroid of the bright core and an outer radius of 
33 pixels also yields P ps 15%. The lack of a signifi- 
cant difference in the polarization between the core of 
the PSF and the much fainter flux surrounding the core 
suggests that the flux in the wings of the PSF completely 
dominates any resolved emission sources in the nuclear 
region. 

The FOC data are not useful for tracking the variations 
in P because no other observations cover this wavelength 
region. What we can conclude from the near-UV imaging 
and polarimetry is that the nucleus and the regions scat- 
tering AGN light are small. We set a conservative upper 
limit for the extent of the scattering regions by using 
the worst of the FOC PSFs, when POL60 was in the 
telescope beam, with FWHM= 3.1 pixel. The data are 
therefore consistent with the near-UV light of the AGN 
(both polarized and unpolarized sources) being contained 
in an unresolved point source with a FWHM— 07044. For 
Mrk 231, this corresponds to a radius of 19 pc. 

Past attempts to resolve the nucleus of Mrk 231 at 
various wavelengths have also yielded only upper limits 
to the size. Ground-based near-IR and 12.5 /mi obser- 
vations, as well as early HST measurements using the 
Near-Infrared Camera and Multi-Object Spectrometer 
(NICMOS) at 1.6 ^m have set limits of -100 pc for 
the FWHM of the unresolved nucleu s (|Lai et alJll99^: 
ISoifer et al.ll2000t loTuiTen et al.ll200l(l . Evidence for a 
rotating molecular torus located 30-100 pc away from 
central continuum source is i nferred by the observa tions 
of hyd roxyl maser em i ssion ( Klo ckner et all 12003). At 
18 cm, lLonsdale et all i|2003|) detected emission extend- 
ing 0"03 southward of an unresolved core component 
with a -4 pc FWHM. Recently, archival NICMOS data 
obtained in 1998 were re-analyzed, taking full advantag e 
of the stable PSF of the instrument ijLow et al J 12004). 
They set an upper limit to the diameter of the core of 
<7 pc at 1.1 [im. As menti oned above , the so phisticated 
PSF analysis employed by lLow et all l)2004t) cannot be 
used for the FOC data because of the large differences in 
the PSFs observed. 

4. RESULTS AND DISCUSSION 

4.1. X-ray Continuum and FeKa Emission 

With the direct X-ray contin uum of Mrk 231 bl ocked 
by a Compton-thick absorber ijBraito et al.l 12004]) . the 
> 3 keV Chandra and XMM-Newton nuclear spectra 
are composites of indirect X-rays that reach the observer 
through reflection off of neutral material and/or scat- 
tering by ionized plasma into the line of sight. In this 
scenario, the FeKa emission is expected to be spectrally 
unresolved at S3 CCD resolution (FWHM-200 eV at 
6.4 keV). Intrinsically broad FeKa emission is generated 
in or near the inner accretion disk; such a line would be 
completely obscured by the same absorber that blocks 
the direct continuum. Our interpretation of the broad 
FeKa feature detected with XMM-Newton by TK as a 



Variability Study of Mrk 231 



7 



composite of narrow lines is more consistent with the 
combined Chandra data (with better photon statistics; 
see Fig. |2J| as well as the physical understanding of the 
system. The best-fitting energies for two narrow lines fit 
separately in the XMM-Newton and combined Chandra 
spectra are 6.3 and 6.6 keV, suggesting reflection off neu- 
tral or nearly neutral material (to generate the lower en- 
ergy fluorescent emission) and scattering off more highly 
ionized plasma (or neutral plasma with an outflow ve- 
locity of ~ 10, OOOkms -1 ). This is generically consistent 
with the preferred model of G02 with multiple scattered 
lines of sight to the nucleus comprising the composite 
2-8 keV spectra. However, at present the relationship 
between the FeKa emission lines and the continuum 
emission is unclear. 

A closer examination of the individual Chandra spectra 
indicates that two narrow Fe Ka lines are not detected in 
all four observations, and therefore these features could 
be flux-variable on timescales of weeks. However, this is 
speculative at present given the photon statistics of the 
individual spectra. 

We confirm the 2-8 keV continuum variability of ~ 
45% within 6 hrs presented in G02 with a reanalysis of 
the Chandra data from observation 1031 using archival 
data and an updated version of the CIAO software. Vari- 
ability is not detected within any of the new observa- 
tions, and the overall 3-8 keV flux variability between 
the 4 Chandra observations is < 15%. The difference 
between the 3-8 keV flux of the XMM-Newton observa- 
tion (2001 Jun) and the combined Chandra observations 
(2003 Feb) is ~+20% over 20 months. This is signifi- 
cantly larger than the measured 2-10 ke V flux differenc e 
between ACIS-S3 and EPIC-pn of - 5% llSnowdenl2 001 1. 
and therefore likely indicates actual long term variability 
in the hard-band flux. 

4.2. Na I D Absorption- Line Variations 

The absorption lin e s invo lving Na I D were organized 
bv iBoksenberg et al.l l)1977l) into redshift systems I (ob- 
served wavelength 6050 A) and II (6015 A), and a system 
at 5984 A was added bv lBoroson et all l|1991fl . Today, the 
latter group is generally known as sys tem III, though it 
was te rmed system IV as recently as iKollatscl mv et 
(1992). 10 All three Na I D systems are apparent in 
the SPOL spectrum of SSAA from 1992-1993 and in 
our 2002—2004 spectrum (both shown at an expanded 
scale in Figure |5J). System Ha remained approximately 
unchanged during the decade spanned by the observa- 
tions, but the system I trough appears to have weakened 
slightly and system III nearly disappeared altogether. 
These effects are not simply a result of a difference in 
spectral resolution, as the spectra were obtained with 
the identical grating (resolution 8 A) , and other features 
in the same spectra, such as the narrow peak atop Ha 
evident in Figure El are as apparent in the 2002—2004 
data. 

The rise and fall of the sy stem III Na I D lines 
have been well chronicled (e.g.. iKollatschnv e ^_aljjl^9_2t 
FRM) . The most recent characterization bv lLipari et al.l 
( 2005) features a smooth decrease in line strength from a 

10 The original system III of |Bokscnbcrg et alj 119771) appears 
only in the Balmer and Ca II H,K lines and is taken to indicate 
the systemic redshift of the host galaxy. 



maximum ( EW = 2.2 A) in 198 8 July to EW = 0.36 A in 
the data of lR.upke et"aT] l|2002fl obtained in 2001 Febru- 
ary. However, data were lacking for a 7 yr interval of 
the decline. Our HET UFOE spectra shown in Figure 
were obtained near the end of the hiatus, yet neither 
Di nor D2 were detected to an individual u pper limit of 
EW ~ 0.03 A, well below the detection of iRupke et alJ 
(2002). At the same time, our spectra from the spec- 
tropolarimeter in 2002—2004 confirm their detection at 
approximately the same strength (Figure |5J. The infer- 
ence is that system III essentially disappeared altogether 
for a period commencing at some point after the obser- 
vations of FRM in 1994 April, was absent during 2000 
January and April, and reappeared weakly by early 2001. 
IRupke et alJ l)2002|) also claim that the observed redshift 
of the feature appears to be increasing with time, though 
much of the historical data do not resolve the doublet. 

Multi-component decomposition of the Na I D com- 
plex has been ca rried out under various assumptions 
l|R,udv et alJlQSfit FRM: lRupke et al.l2002l) . and will not 
be repeated here. We simply note that, apart from the 
absence of system III, th e HET spectra closely resemble 
the data of IRupke et alJ (|2002fl obtained one year later. 
The higher resolution of the UFOE data allows a closer 
examination of the narrow features near the edges of the 
principal system I feature. The lie doublet, whose Di 
component lies at the blue edge of the main system I 
trough, has progressively increased in strength over the 
past decade. 11 In fact, the now-obvious D2 component 
lies in what was a region of the continuum in the 1991 
data of FRM and was only beginning to appear in their 
1994 spectrum. There is marginal evidence that system 
lie D2 is slightly narrower than the system Ha features, 
as indicated in Tabled however it is clear that the as- 
signment of lines to distinct systems is meaningful only 
for bookkeeping purposes. There is also marginal evi- 
dence that system la Di and Ii D2 are each composed of 
two narrow components with a separation of ~1 A. We 
consider the identification of the Ii Di component at the 
blue edge of the main system I trough to be tentative, 
and a higher signal-to-noise ratio spectrum would be use- 
ful for separating components in this region. Finally, we 
point out a curious excess in the "continuum" immedi- 
ately surrounding the system Ha Di line in Figured A 
similar feature is apparent at somewhat lower signal-to- 
noise ratio in the 1994 spectrum of FRM; together the 
observations suggest the existence of a narrow emission 
component of Na I D associate d with one of the a bsorp- 
tion systems (see also FRM and lLipari et al.ll994l on this 
issue) . 

As pointed out by FRM, an interpretation of the main 
syste m I trough in term s of a few highly saturated dou- 
blets (jRudv et alJ T9851 is not feasible because a column 
density in neutral Na as high as 10 15 cm -2 develops 
strong damping wings, inconsistent with the observed 
profile. FRM's choice of a number of optically thin com- 
ponents (An & ~ 10 13 cm" 2 ) over several moderately sat- 
urated doublets is based on the presence of residual flux 
at the bottom of the trough. However, they appear to 
have ignored the underlying galaxy starlight, which our 
imaging has shown is significant in the typical observing 
aperture of a ground-based telescope. The general flat- 

11 System lib has not been seen since the 1991 spectrum of FRM. 
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ness of the system I trough suggests that the individual 
components are optically thick, so it would seem that a 
blend of a number of moderately saturated lines is still an 
option, and the individual system column densities may 
well be intermediate between the above two estimates. 

Starlight as the sole source of residual light in the sys- 
tem I trough is consistent with the lack of polarization 
in the bottom of the trough in 1992-1993 (Figure EJ, 
and SSAA estimated the galactic dilution of the AGN 
polarization accordingly. From Figure we also find 
that the residual flux level in the echelle spectrum of the 
trough, ^17% of the surrounding continuum level, is con- 
sistent with the host galaxy dilution expected within a 2" 
aperture. In contrast, the 2002—2004 spectropolarimetry 
indicates measurable polarization in the trough, with a 
position angle common to that of the surrounding contin- 
uum. To be sure, the net polarization in this feature is a 
factor of two below that of the surrounding continuum, 
even considering the decreased continuum polarization 
measured in 2002—2004. However, its presence suggests 
that some scattered light may have been introduced into 
the trough by 2002. A small amount of scattered AGN 
light that avoids absorption in the main system I trough 
would tend to support our tentative identification of the 
Ii Di feature above. 

A quantitative estimate of the intrinsic degree of po- 
larization of the scattered light in the system I trough 
of Mrk 231 is strongly affected by the uncertainty in the 
host-galaxy contribution within the observing aperture, 
but it could exceed the polarization of the surrounding 
continuum. In high-z BAL quasars, scattered light com- 
pletely dominates the residual flux in the rest-frame UV 
troughs, and the measured polarization c an exceed 10% 
HScnrnidt fe HinesilT999t l Ogle et al .1 179991). 

4.3. Nature of the Optical Polarization Change 

The magnitude of the polarization observed in 
Mrk 231 and the similarity of characteristics between 
the emission line and continuum polarization have long 
been interpreted as indicating an origin in s cattering 
l|Thompson et al.lll980t ISchmidt fe Millerlll985D . In red, 
strongly polarized quasars in general, and Mrk 231 in 
particular, there is compelling evidence in the shapes 
of the polarized and total flux spectra that the scat- 
tering occurs in exten sive dust clouds located within 
~100 pc of the nucleus llGoodrich fe Milledll99i SSAA; 
Mines et all l200lt iSmith et alJ 120031 l2004|) Again spe- 
cific to Mrk 231, our HST imaging polarimetry pre- 
sented above places these clouds interior to a 20 pc ra- 
dius, and the fact that polarization position angle struc- 
ture is observed across the Ha emission profile f H3.2l 
SSAA) indicates that the broad-line region is partially re- 
solved from the vantage point of the scattering medium, 
i.e., much of the scattering could occur far deeper in 
the nucleus. Th e 1.1 /^m nuclear light constrained by 
lLow et all l)2004|) to originate inside a radius of ~3 pc 
is apparently emitted by hot dust grains, the survival 
of which is possible as close as <1 pc to a source with 
L ~ 4 x 10 12 Lq. In view of these facts, we imagine 
that all of the optical emission escaping from the nucleus 
may have undergone at least one scattering. Indeed, 
with optical extinction in interstellar clouds scaling as 
ATh [A y ~2x 10 21 cm' 2 mag' 1 (e.g.. lBurstein fe Heilesl 
1978), the very high absorption measured toward the X- 



ray power law (N H - 2 x 10 24 cm~ 2 ; Bra ito et alJ ktMl 
would suggest that there are no open sightlines to the 
central source from soft-X-r ay to submillimeter w ave- 
lengths. However, as noted by Maioli no et alJ 1)200 lj) . the 
column densities from X-ray and optical extinction mea- 
sures in AGN often imply much lower dust-to-gas ratios 
than typical interstellar clouds. Furthermore in Mrk 231, 
the X-ray absorber may be largely dust-free; it is likely to 
be much more compact than both the longer wavelength 
continuum emission and the dusty shroud obscuring the 
optical-UV continuum source. 

Because scattering in a clumpy dust shroud is prone 
to the searchlight effect (redirection of the illuminating 
beam), it is not possible to relate directly the time scale 
for a polarization change to motion of clouds. However, 
the orbital period for gas at a distance of 0.5 pc from a 
10 9 Mq black hole is 1000 yr, and so a significant change 
in scattering geometry over a decade is not implausible. 

The following essential phenomena must be explained 
with regard to the observed polarization change in Mrk 
231: 1) a decline in degree of polarization P, and in 
polarized flux P x F\ by a factor of ~l/3 (e.g., P bs = 
3% to 2% at A - 7000 A) in both the continuum and 
the emission lines; 2) a position angle of polarization 
that has considerable structure across the spectrum but 
is manifestly constant with time; 3) an upper limit on 
the amount of any systematic change in total nuclear 
brightness of 15%; 4) a spectrum of polarized flux that 
has narrower emission lines than that in total flux; and 
5) a reduction in the depth of the system I trough as well 
as the appearance in 2002 of significant polarization in 
this feature at the systemic polarization position angle. 

4.3.1. The Importance of Starlight Dilution 

An increase in the relative dilution by an unpolarized 
component is an obvious means of decreasing both the 
degree of polarization, P, and the polarized flux, P x Fx, 
without modifying the position angle. As an example, a 
decline in the intrinsic brightness of the AGN, while the 
diluting host galaxy starlight remains constant, would 
lead to a drop in the polarization measured within a fi- 
nite aperture. The allowed magnitude of the decline is, 
of course, subject to our conclusion that the overall nu- 
clear brightness varied by no more than 15% during the 
decade. 

The lack of an aperture dependence to the polariza- 
tion measurements summarized in Table|3]aheady argues 
that starlight cannot be a strong diluting component to 
the light measured through a typical spectropolarimet- 
ric aperture and thus is not an important factor in the 
polarization change. The strongest evidence for this con- 
clusion is provided by the residual flux in the base of the 
system I Na I D trough, measured at 15% in 1993 and 
21% in 2003 in the 3" aperture used with the spectropo- 
larimeter. With stars initially accounting for no more of 
15% of the total light at 6000 A, a change of nearly 1/3 
in polarized flux at that wavelength requires the AGN to 
have faded by more than a factor of four. The net spec- 
trum in 2002 would have been far fainter than what was 
measured and the resulting Na I D trough would have 
been <60% of the continuum level. 

Additional evidence of a change in polarization of the 
AGN component itself comes from the fact that the in- 
trinsic (AGN-only) _B-band polarization determined from 
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the imaging polarimetry in 1996 (Pb = 6.45%, Figure lTUjl 
is no larger than the net (AGN + galaxy) value computed 
from the spectropolarimetry for the same filter passband 
in 1992—1993 (6.62%). Indeed, the imaging polarime- 
try implies that the S-band polarization appropriate to 
a 3"xl0" aperture in 1996, 5.70% (Table EI), is already 
well on its way to the 2002-2004 value of 5.04%. 

Finally, we note that increased dilution by the opening 
of a new, direct (i.e., unscattered) sightline to the AGN 
could avoid the problem with the depth of the Na ID 
trough if the direct beam also passed through the sys- 
tem I clouds. However, this mechanism would leave the 
polarized flux unchanged, contrary to what is observed. 

4.3.2. A Strongly Polarized Component 

An important clue toward a more successful model for 
the polarization change was noted in M3.2I namely, the 
spectrum of polarized flux is not simply a reflection of the 
total flux spectrum. Specifically, the individual Fe II and 
H I emission features are much more sharply defined in 
the former. This is not simply a matter of contrast, but 
more of line width, and implies that the scattering clouds 
do not see the full velocity dispersion of the broad line 
region. The physical origin of the difference will be ad- 
dressed below, but here we point out that the difference 
between the two spectra requires that a large fraction of 
the AGN total flux spectrum, containing the continuum, 
Na I D absorption, and smeared emission lines, has no 
polarization whatsoever. This leads to an upper limit on 
the amount of light that becomes polarized by scattering 
(because scattering is not 100% efficient at inducing po- 
larization), and therefore to a lower limit on the degree 
of polarization of the scattered light. 

We have carried out simple experiments to estimate the 
fraction of the total flux spectrum from 2002—2004 that 
is scattered by simply summing the spectrum of P x F\ 
and F\ in various relative fractions and judging by eye 
when the sharper features of the former are no longer 
visible in the total. Conservatively, we find that a higher 
degree of structure would be apparent in the total flux 
spectrum if a spectrum like that of the polarized flux 
comprised more than 30% of the total at a wavelength of 
5500 A. The implications are that the more sharply struc- 
tured, polarized AGN component must be diluted by a 
virtually unpolarized AGN component (and starlight) to 
a ratio of at least 2.3:1, and that the intrinsic degree of 
polarization of the scattered light at this wavelength in 
2002 was a substantial P sca > 3.3 x P ohs = 12.5%! This is 
comparable to the observed polarization level at blue op- 
tical and near-UV wavelengths (SSAA), suggesting that 
the unpolarized AGN component is primarily responsible 
for the overall very reddened optical-UV spectral energy 
distribution, and that its diluting effects are greatly re- 
duced at short wavelengths because of dust extinction. 
Because the polarized component is such a small frac- 
tion of the total light at visible wavelengths, a drop in 
the net observed polarization from 4.8% to 3.8% between 
the 1990's and 2000's at 5500A only requires a decline 
in brightness of the polarized AGN component of <30%. 
The change in brightness of the overall nucleus (AGN + 
starlight) is then <14%, consistent with the limit on the 
brightness change between the epochs. 



4.3.3. More Complex Mechanisms 

Alternate explanations for the polarization change in 
Mrk 231 can be constructed by the introduction of ad- 
ditional polarized components between the epochs. We 
have investigated this option by attempting to synthe- 
size the observed flux and polarization spectra for the 
two ep ochs using the method applied bv I Schmidt et all 
(2002) to the bipolar protoplanetary nebula surround- 
ing CIT 6. We find that such a model can reproduce 
a large amount of the polarization structure seen in the 
Mrk 231 spectrum, and that the observed changes can 
be interpreted as a change in the relative brightness of 
the polarized components. For example, the opening of 
a new scattered sightline with a polarized flux spectrum 
like the difference spectrum shown in Figure |H1 but with 
a polarization position angle orthogonal to the systemic 
angle, i.e., 9 ~ 10°, provides an excellent match to the 
most recent observations. If the fractional polarization 
of the new sightline were large, P > 10% at 6500 A, the 
overall increase in brightness through our spectroscopic 
aperture would be <15%. At some level, changes in 9 be- 
tween epochs are also expecte d from this interpreta tion. 
Unfortunately, as discussed bv lSchmidt et all (|2002|) and 
SSAA, even simple two-component models suffer from 
being non-unique, with a range of combinations of po- 
larization, position angle, and flux level reproducing the 
observations with equal quality. 

4.4. Dust- Reddened Quasars and the Structure of AGN 

Mrk 231 can be considered the first known example 
of a population of optically red, IR-luminous AGN with 
quasar-like total energy outputs that are now being iden- 
tified in significa nt numbers by near-IR, optical and X- 
ray surveys (e.g.. Mines et al.ll200H iSmith et al1 f2002b; 
!Mattl2002i iZakamska et al.l2005l and references therein) . 
These objects are found in both Type 1 (broad permit- 
ted lines) and Type 2 (narrow permitted lines) varieties, 
and are distinguished from other samples of quasars by 
the fact that a sizeable fraction of each type (>10%) 
show optical polarization of the total light P > 3%, the 
traditional threshold for "highly" polarized AGN. Gen- 
erally, the polarized Type 2 quasars have been shown 
to follow the basic predic tions of the Unified Scheme 
for AGN l)Antonuccil 119931) by revealing broad permit- 
ted lines characteristic of a Type 1 object in polarized 
light. This indicates the existence of a broad line region 
and attests to anisotropic obscuration surrounding the 
centr al direct continuum source (e.g., IZakamsk a et alJ 
12005(1 . The optical obscuring structure is nominally re- 
garded as a dusty /m olecular torus at r ss 1-100 pc (e.g., 
iNenkova et alJl2002tl . which is seen at high inclination 
(i.e., edge-on) in the Type 2 systems, and the polariza- 
tion is produced by scattering off particles located near 
the polar axis. 

Though the Unified Scheme has been very effective in 
accounting for many key aspects of the various empirical 
classes of Seyfert-luminosity AGN within a single simple 
framework, it provides no obvious explanation for the 
strongly polarized Type 1 Seyfert galaxies and quasars 
like Mrk 231, which by virtue of the strong broad- line 
component in their total flux spectra are thought to be 
viewed at smaller inclination angles — i.e., through the 
"hole" of the torus. In total flux, the new population of 
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red quasars shows evidence for considerable reddening, 
whereas the p olarized flux spectra are considerably bluer 
iSmith et alJl2003fl . 

In an attempt to account for the ~1% polar i zation 
often seen in Typ e 1 AGNs. iGoodrich & Milled l|1994fl 
and more recently iSmith et alJ l)2002al 120051) proposed 
the existence of an equatorial disk nominally coplanar 
with and lying inside the torus. Scattering off dust 
and/or electrons in this disk provides an additional po- 
larizing mechanism for Type 1 objects. Depending on 
inclination angle, the observed polarization of a Type 1 
AGN might be dominated by scattering off particles 
near the polar axis (the predominant Type 2 polarizing 
mechanism, for i ~ 45°) or off the e quato rial disk (for 
i < 45°). The scheme of ISmith etafl <|2005j) is currently 
a geometrical framework rather than a physical model, 
but there are clear parallels with disk-lau nched winds 
(IKonigl fc Kartiell994UMurrav et alilQQl . If the equa- 
torial scattering medium is identified with the disk wind, 
then this mechanism should not be an important source 
of polarization in BAL quasars such as Mrk 231. 

4.5. Toward a Model for the Nucleus of Mrk 231 

Based on the multiwavelength data available on 
Mrk 231, the nucleus is clearly complex with more than 
one structure scattering continuum and emission-line 
photons into the line of sight. Insight into the nucleus 
on the smallest physical scales in Mrk 231 is provided by 
X-ray observations. The most interesting constraints are 
provided by the combination of X-ray data from differ- 
ent observatories. The confirmation of Compton-thick 
absorption from t he detection of the direct continuum 
at hard X-rays bv lBraito et al.l l)2004|) coupled with the 
continuum variability detected by G02 point towards a 
compact Compton-thick absorber. This opaque absorber 
must lie interior to the broad-line region, and therefore 
would not be identified with a molecular torus as are the 
Compton-thick absorbers in Type 2 AGNs. The X-ray 
absorber is likely to be on scales smaller than or com- 
parable to the UV-optical continuum region, to enable 
detection of the continuum and radiative driving of the 
BAL outflow. As discussed in G02, the variability of 
indirect 2-8 keV X-rays constrains the scale of the scat- 
terer to be ~ 10 15 cm -2 and consequently the ionization 
of this plasma is expected to be high (i.e., with Fe al- 
most totally stripped of electrons). The new combined 
Chandra data indicate that highly ionized plasma may 
have been detected in emission (with one Fe Ka emission 
component at ~ 6.6 keV), though the relationship be- 
tween the Fe Ken-emitting gas and the observed 2-8 keV 
continuum emission is unclear. The dominant scattering 
mechanism in this case is electron scattering, and so this 
structure is distinct from the dust-enshrouded scatterer 
dominating the optical polarization. 

On larger scales, two polarization characteristics sug- 
gest that there exists a prevailing symmetry in the un- 
resolved active nucleus, and that some of the structures 
inferred by the Unified Scheme may be relevant. First, 
the constancy of the polarization position angle over 3 
decades, despite a sizable change in P, requires that the 
dominant optical scattering region(s) maintain the same 
orientation relative to the continuum source and broad 
line region. Second, the individual Fe II and H I emis- 
sion features are much more sharply defined in the po- 



larized flux spectrum than in the total flux spectrum. 
The difference is better displayed in Figure ITT1 where we 
compare the total flux (less the contribution of the host 
galaxy; thick line) and polarized flux spectra (thin lines) 
for the observed- frame 5000— 6000 A region, with various 
amounts of gaussian smoothing applied to the polarized 
flux. As can be seen, the spectrum of P x F\ must be con- 
volved with a gaussian of FWHM = 1500-2000 km s" 1 
to approximate the smoothness of emission features in 
the total flux spectrum. Also apparent in Figure 1111 is 
the lack of perfect correspondence between the velocity 
centroid of features in the two spectra, with discrepan- 
cies as large as ~500kms _1 evident between the two. 
Some of these displacements are likely associated with 
the small polarization position angle rotations observed 
across the profiles of prominent emission features. 

An attractive explanation for the narrower spectral 
features in the polarized flux spectrum is that the scat- 
tering plane is essentially orthogonal to the direction of 
an important dynamical motion of the line-emitting gas. 
Scattering off particles in polar lobes, as in the Unified 
Scheme for AGNs, offers such an arrangement, if the 
broad-line gas orbits the nucleus in a disk whose axis is 
nominally coincident with the lobes. The lobes also pro- 
vide a natural explanation for the very large polarization 
inferred for this component in H3.2I The broad line re- 
gion disk would also be roughly coplanar with the dust 
torus that presumably determines whether we perceive 
a Type 1 or 2 object in the total flux spectrum. Since 
Mrk 231 is a Type 1 AGN, the dust torus is not fully 
edge on, but it still serves to define the polar scattering 
region(s) at a position angle of ~7°. The polarization 
properties are also consistent with an equatorial BAL 
outflow, though the covering fraction of the BAL wind is 
apparently greater than that of the torus. Furthermore, 
the detection of some polarized continuum unabsorbed 
by Na ID (FigJZJl shows that the wind does not cover 
the entire sky. 

The polarization position angle is closely aligned with 
the axis of the 40-pc sc ale triple source (core- lobe) im- 
aged at GHz frequencies (|Ulvestad et al.ll999al and refer- 
ences therein). On parsec scales, a double radio source in 
the n ucleus is misaligned by ~ 65° from the larger struc- 
ture l|Ulvestad et al.ll999bl) . and therefore not coincident 
with the position angle of the scattering medium. Ignor- 
ing any velocity dispersion within the scattering lobes, 
differential projection effects across the lobes, and the 
unknown inclination of the orbiting disk to the line of 
sight, a velocity of 1500—2000 km s _1 corresponds to 
gas orbiting 1—2 pc away for a black hole at the high 
end of the mass distribution, 10 9 Mq. Even this upper 
estimate for the radius is we ll within the upp er limit to 
the AGN size determined bv lLow et"aTI l)2004|) . A broad 
line region with dimensions of a few pc is larger than 
the ^0.3 pc estimated for objects of comparable lumi- 
nosity (XL ° ~ 10 45 ergs s _1 ) using the reverberation 

technique : Kaspi et alJl2000(l . but the assumptions made 
above suggest that the velocity dispersion estimate is a 
l ower limit, and so th e disk size may be overestimated. 

ISmith et al.l l)2004|) also consider the specific case of 
Mrk 231 and, based on the redness of the optical con- 
tinuum, the presence of BAL features, and the general 
constancy of the polarization position angle across the 
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spectrum, conclude that scattering in polar lobes is the 
dominant polarizing mechanism. However, their assign- 
ment to "far-field" scattering in the lobes is not sup- 
ported by our much higher signal-to-noise ratio data that 
reveal considerable position angle structure across the 
Balmer and Fe II emission lines, nor by the HST po- 
larimetry result that locates the scattering material no 
further than 2 pc from t he con tinuum source. In fact, 
the results of ILow et all l|2004|) suggest that the scat- 
tering in the polar lobes likely takes place within just a 
few pc of the nucleus; a much more compact scattering 
region than typically assumed in the Unified Scheme or 
actually observed in se veral nearby Seyfert g alaxies and 
some quasars (see e.g..fZakamsk a~et al.ll2005l) that show 
scattered light from regions of the order >1 kpc. How- 
ever, light polarized by scattering from such a compact 
region would likely also be blocked by the same structure 
that blocks the broad-line region (with similar scales) in 
Type 2 objects. 

The observations also require that the compact pri- 
mary scattering region be enshrouded in dust to account 
for the lack of a narrow emission-line region. In this view 
of Mrk 231, the much reduced size of its polar scattering 
regions and the fact that light from these regions must 
also traverse a significant amount of dust are significant 
departures from the simple unified picture of AGNs. The 
compactness of the scattering regions and the ubiquitous 
dust distribution around the nucleus and possibly within 
the polar lobes may, in fact, be responsible for the large 
decrease in P observed between the early 1990s and early 
2000s. One possible scenario is that dusty material in- 
terior to the main scattering clouds changed the amount 
of illumination within the scattering lobes or, alterna- 
tively, dust exterior to the polar lobes blocked more of 
the scattered light along our line of sight. In either case, 
large flux changes are not necessary to cause the polar- 
ization variation, and the position angle of the polarized 
flux will be prese rved. Finally , the a ugmented morphol- 
ogy proposed bv iSmith et al.1 l)2004D . i.e., an additional 
scattering disk coplanar with the dust torus, but interior 
to it, may be responsible for the variation in polarization 
position angle across emission features and the overall 
rotation of 9 observed throughout the optical and UV 
(SSAA). 

4.6. Future Investigations 

The holistic examination of the indirect X-ray, UV, 
and optical emission has provided some new insights into 
the nuclear structure of Mrk 231 that bear further inves- 
tigation. From the close examination of the individual 
Chandra spectra, it appears that two narrow Fe Ka lines 
are not equally strong in all four observations. A claim 
of flux variability between the Fe Ka lines is speculative 
at present given the photon statistics of the individual 
spectra. With its higher effective area at E > 3 keV 
compared to Chandra, additional ~ 30 ks XMM-Newton 
observations with similar time sampling to these Chan- 
dra observations would be sensitive to this phenomenon. 



A detection of Fe Ka line variability on these time scales 
would constrain the location of scattering material to 
within the broad line region, as suggested by the con- 
tinuum variability detected by G02 (and assuming the 
scattering media are also the sources of the Fe Ka emis- 
sion). A single, long-look (100-300ks) XMM-Newton 
observation would also be productive for searching for 
additional hard-band variability on intermediate time 
scales of —day, for studying the iron-band features in 
depth, and for potentially detecting soft-band emission 
lines from the starburst and AGN with the Reflection 
Grating Spectrometer. 

The lack of [O III] emission in the optical spectrum 
indicates that the covering fraction of dust around the 
nucleus is very high, thereby preventing ionizing pho- 
tons from illuminating low-density gas in what would 
otherwise be the narrow-line region. The large-column 
density X-ray absorber clearly does not cover the entire 
sky, otherwise short-timescale variability would not be 
detectable. As X-rays can easily penetrate dusty gas that 
would block FUV photons, the X-ray narrow-line region 
(ionized by X-ray-energy photons) should exist and may 
even have already been detected (TK). In this region, 
the strongest features expected are the forbidden lines o f 
O VIIA22.1A and NeIXA13.7A fe.g.. lQgle etafl fcOO.U 
High-resolution X-ray spectroscopy with the next gener- 
ation of X-ray observatories such as Constellation- X will 
enable a detailed emission-line study undiluted by the 
strong direct continuum, a rare opportunity to study the 
X-ray narrow-line region of a nearby, quasar-luminosity 
AGN. 

The polarization variation we have detected in Mrk 231 
is perhaps the most accurately recorded so far among 
radio-quiet AGN, but it is not unique. I Zw 1, which 
also exhibits a polarized spectrum caused by scattering, 
displayed a reduction in fractional amount and a rotation 
of —50° in po sition a ngle between the 19 80's and 1990's 
( Smit h et alJll997|) . iGoodrich fe Millerl ()1994j) noted a 
25% decline in the polarization of Mrk 509 over a 2 yr 
period, an d a similar va r iation i n the same object wa s dis- 
cussed bv lYoung et alJ l|1999lh ISmith et all l|2002aD also 
point out apparently significant changes in Akn 120 and 
NGC 7469 over periods of -18 months. The fact that all 
are Type 1 AGNs suggests that a significant fraction of 
the polarization of these objects arises within the inner 
few pc of the nucleus, and offers hope that further stud- 
ies of polarization and its variability might continue to 
improve our understanding of their nuclear structure. 
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TABLE 1 

Chandra Observing Log 



Obs. ID 


Date 


Exposure 


Nuclear Count Rates 


a (10~ 2 ct s- 1 ) 






(ks) 


0.5-2.0 keV 


2.0-8.0 keV 


1031 


2000 Oct 19 


39.3 


1.98 ±0.07 


1.92 ±0.07 


4028 


2003 Feb 03 


39.7 


1.86 ±0.07 


1.98 ±0.07 


4029 


2003 Feb 11 


38.6 


1.83 ±0.07 


2.00 ±0.07 


4030 


2003 Feb 20 


36.0 


1.79 ± 0.07 


2.04 ±0.08 



a The values are background-subtracted, though the background was negligible in the 
hard band and < 3% in the soft band. Nuclear source counts were extracted from a 
3"-radius source cell. Background counts were collected from an annulus with inner and 
outer radii of 5" and 10", respectively, and normalized to the area of the source cell. 



TABLE 2 
Chandra X-ray Spectral Fitting 21 











Observation Identification 






abi 


Parameter (Units) 


1031 


4028 


4029 


4030 


4028-30 


XMM-Newton h 


XMM-Newton c 


VI 


3-8 keV Count Rate (10~ 2 s" 1 ) 


1.32±0.06 


1.38±0.06 


1.52±0.06 


1.52±0.07 


1.51±0.04 






CO 


5.8-6.7 keV Count Rate (10~ 2 s" 1 ) 


0.21±0.02 


0.18±0.02 


0.23±0.03 


0.19±0.02 


0.19±0.01 






pi 


F 3 _ 8 (10" 13 crgcm- 2 s- 1 ) 


4 fi5+ - 29 


4.80l°;g 


c 17 +0.35 
- i ' -0.35 


5 05+ 37 
°- uo -0.37 


5 03+ ' 20 


4 lfi+0- 26 
4 ' iD -0.27 


4 16+ ' 27 

4 - 1D -0.26 




L 3 _ 8 (lO^ergs- 1 ) 


i OO+0.09 
l-83_ 14 


i-87±8:i| 




i-98±°:ii 


i 07+O.O8 
L - v ' -0.09 


1 60+° 10 


i cc+0.10 

1.66I n 


c 


r 


o.47i8:g 


o.48i8:2? 


•W+ ' 40 


o 51+0- 24 

U - O1 -0.49 


n41 +0.20 
u - 4i -0.24 


0.83l8;l 


o.76l8; 3 2 8 




Power Law Normalization (10 -5 /E d ) 


2 .2oit;g 


2-ao±i:g5 


2 0Q+ 1 ' 63 


r) qq + 1.22 


2 31 +0.79 


o 71 +2.56 
°- ' 1 -2.20 


3.45+?6? 


!*r 


Line Energy (keV) 


fi n 4 +>0.66 


6.68l°;8? 


6 ^4+ 04 
"■^-0.04 


6.4 


6.311H? /6°64lH 5 7 


6 47+O- 13 
D - 4 '-0.13 


6.3il8;ii/6.62l8;i? 


to 
cc 


cr (keV) 


45+ ' 72 


0.01 


0.01 


0.01 


0.01/0.01 


19+° 12 
u - ±M -o.io 


0.01/0.01 


i— ' 


Line Flux (10 -6 photcm~ 2 s _1 ) 


6.00+ 4 6 7 ° 


2.8ilr° 7 ? 


3.091 2 ;" 


< 1.86 


n qq+0-73,, 2C .+0.97 
u - y,;> -0.84/ i - zo -0.81 


3.56l?;|l 


1.65li;8 2 9 /1.73li; 4 




Line EW (eV) 


638li™° 


28ll?? 4 




< 160 


70lg/98l 77 


483l 2 ? 3 


143lf 5 /167li8 8 5 




Edge Energy (keV) 


7.12 


7.12 


7.12 


7.12 


7.12 


7.12 


7.12 




^raax 


29+ 71 


ii+ - 76 


49+0-70 


17+0- 75 

u - i ' -0.17 


o^l^ 


2 o+o.3i 

u '^°-0.23 


0.2818;^ 




x> 


23.2/27 


22.4/29 


21.7/31 


13.1/31 


60.1/72 


23.0/27 


21.6/26 






0.68 


0.80 


0.89 


0.99 


0.65 


0.67 


0.71 





a Errors are for 90% confidence taking 1 parameter to be of interest (Ax 2 = 2.71). Values with no errors were fixed for the fitting. The errors for -Fb—g, the observed-frame 3—8 keV flux, 
and L 3 -8, the rest-frame 3—8 keV luminosity, were determined from 90% confidence errors in the power law normalization with all other parameters frozen to the best-fit values. 
b Thc XMM-Newton EPIC-pn data from 3.0—10.0 keV were fit for this analysis. These spectral-fitting results are from the same data set presented in TK, and are fully consistent. 
c The same data, fit with two unresolved Gaussians instead of a single Gaussian with a free to vary. 
d Units of /e are phot cm -2 s -1 keV -1 . 
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TABLE 3 

Log of Optical and UV Observations 



Date 


Exposure 


Coverage 


AA 


Aperture 


Telescope 


pa 


6» a 


(yyyymmdd) 


00 


(A) 


(A) 


(") 




(%) 


(°) 


Spectropolarimetry 


First Epoch (1990's) 
















19920312 


800 


4200- 


-7900 


17 


3x 10 


Bok 


3.36 


97.5 


19930226 


1200 


4100- 


-7400 


20 


4x 8 


Bok 


3.41 


96.9 


19930428 


1600 


4000- 


-7300 


17 


3x 8 


Bok 


3.40 


96.6 


19930513 


4800 


5630- 


-7360 


8 


3x 9 


Bok 


3.46 


96.6 


Mean values. . 














3.41 (0.04) 


96.9 (0.4) 


Second Epoch 


(2000's) 
















20021204 


2880 


4200- 


-8200 


17 


3x 10 


Bok 


2.45 


96.7 


20030202 


4800 


4200- 


-8200 


17 


3x 7 


Bok 


2.65 


98.1 


20030203 


2400 


4200- 


-8200 


20 


4x 14 


Bok 


2.42 


97.3 


20030203 


12800 


5250- 


-7700 


8 


3x 13 


Bok 


2.33 


97.0 


20030307 


3200 


4200- 


-8200 


17 


3x 10 


Bok 


2.57 


97.3 


20040426 


960 


4200- 


-8400 


16 


0.7x 5 


MMT 


2.55 


98.3 


20040426 


960 


4200- 


-8400 


19 


l.lx 5 


MMT 


2.55 


98.2 


20040426 


960 


4200- 


-8400 


22 


1.4x 5 


MMT 


2.56 


98.4 


20040426 


960 


4200- 


-8400 


24 


1.8x 6 


MMT 


2.56 


98.8 


Mean values. . 














2.52 (0.10) 


97.8 (0.7) 


Imaging Polarimetry 


19961207 


960 


B 


-1000 


3x 10 


Bok 


5.70 


93.4 


19961207 


320 


R 


~1600 


3x 10 


Bok 


2.60 


95.9 


19981128 


15520 


F346M 


~430 


0.47 


HST 


15.3 


97.8 


Echellc Spectroscopy 


20000106 


3x200 


4534- 


-9154 


0.8 


2 (fiber) 


HET 






20000113 


4x900 


4534- 


-9154 


0.8 


2 (fiber) 


HET 






20000403 


3x900 


4416- 


-8810 


0.8 


2 (fiber) 


HET 







a For the spectropolarimctric observations, the mean polarization properties are com- 
puted over the observed range 6100— 6600 A, a largely continuum band common to all 
of the datasets. For a single dataset, uncertainties are dominated by calibration, see- 
ing, and other aperture-dependent effect s as indicated by the standard deviations (in 
parentheses; AP ~ 0.1%, A8 ~ 1°). See H3.2l for further discussion. 



TABLE 4 

NA I D ABSORPTION IN MRK 231 



System 3, 


Ao 


A b 


V 


V - V sys C 


FWHM 


EW 




(A) 


(A) 


(kms" 1 ) 


( kms ) 


( kms -1 ) 


(A) 


I 


5889.950+5895.923 


6047 


7750 


-4650 


~ 1100 


~ 18 


la 


5895.923 


6058.82 ± 0.05 


8168 ± 3 


-4232 


110 ±8 


0.63 ± 0.12 


lb 


5889.950 


6037.97 ± 0.07 


7440 ± 3 


-4960 


101 ± 8 


0.34 ±0.11 


Ha 


5889.950 


6013.49 ± 0.03 


6222 ± 2 


-6178 


92 ± 6 


0.54 ± 0.02 


Ha 


5895.923 


6019.60 ± 0.03 


6223 ± 2 


-6177 


81 ±8 


0.34 ± 0.02 


lie 


5889.950 


6026.81 ± 0.10 


6885 ± 5 


-5515 


73 ±9 


0.28 ±0.03 


lie 


5895.923 


6033.28 ± 0.12 


6903 ± 6 


-5497 


48 ± 15 


0.13 ± 0.04 



a The system identifications are adopted from FRM. 

b Uncertainties were estimated from Monte Carlo simulations using the noise measured 
in the continuum. 

c The measured radial velocity of the line system relative to the systemic redshift of 
12,400 kms" 1 . 
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Fig. 1. — Count rates for the four different Chandra observations of Mrk 231. The decrease in the 0.5-2.0 keV count rate from observation 
1031 to 4028 is likely a result of the declining soft-band effective area over time. Though the full-band count rate has remained remarkably 
stable, the slight changes in soft and hard-band count rates suggest a hardening of the spectrum. 




Fig. 2. — X-ray spectrum and model from the combined 3-8 keV Chandra data from observations 4028—30. The thick histogram is the 
best-fitting power-law model with two narrow gaussian emission lines and an Fe edge (see Table convolved with the Chandra response. 
Vertical dashed lines mark the best-fitting energies of two narrow Fe emission lines. 
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- Residual plots for the four 3-8 keV Mrk 231 spectra from joint-spectral fitting. For the joint-fitting, the values of F and T max 
power law normalizations were free to vary, and the edge energy was fixed to 7.12 keV, the energy of the neutral Fe K edge. To 
the continuum fit, the iron-band region was ignored. The energies of the two possible Fe lines arc marked with vertical dotted 
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Fig. 4. — Contour plots of the line flux versus rest-frame energy of the iron lines detected in the two middle observations. The contours 
are 68% and 90% confidence for two parameters of interest (Ax 2 = 2.31 and 4.61, respectively). Though the fluxes and EWs are consistent, 
the energies arc not. 
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Fig. 5. — Left panel: Full-band, csmoothed image of Mrk 231 from the merged data sets of the three new observations. North is up and 
east is to the left. This 114 ks exposure reveals ring structure around the nucleus as well as a new point source to the southwest of the 
nucleus. Right panel: Three-color csmoothed image of Mrk 231 from the merged data sets of the three new observations. Red, green, and 
blue correspond to 0.35—1.2, 1.2—3.0, and 3.0-8.0 keV, respectively. The blue emission is slightly asymmetric, with an enhancement just to 
the west of the nucleus. The flat nuclear spectrum leads to the white appearance of the nucleus. 
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Fig. 6. — Radial surface brightness profile in the 0.35—2.0 keV b and of the nucle us of Mrk 231 from the merged data of observations 
4028-30. The la errors have been d etermined using the method of IGehrela $1986). The vertical, dotted lines indicate the radii of the 
optical shells, SI and S2, discussed by Li'pari et al. 12005). 
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Fig. 7. — The Na I D absorption systems of Mrk 231 during 2000. The averag e spect r um o f 2000 January and April is shown with 
individual line systems labeled. Systems I and II were first identified by IBoksenberg et alj 119771) . and are seen to break into subsystems 
la, Ii, and He in these much higher resolution HET echelle observations (see also FRM). The Di component of system Ii is tentatively 
identified within the main trough of system I at about the correct position relative to D2. Positions of the system III lines are marked, 
but they are not detected in these observations of early 2000. Shown below the average spectrum is the difference spectrum (Jan — Apr), 
which shows little evidence for change in the ~4 month interval. 
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Fig. 8. — Comparison of the observed-frame, co-added polarization spectra of Mrk 231 from 2002—2004 (thick line) and 1992—1993 (thin 
line). Note the systematic decline by nearly 30% (AP = 1 — 2% across the bandpass) in the degree of polarization with largely unchanged 
position angles. The polarized flux panel (P X F\) also shows the difference spectrum between the two epochs (1990's minus 2000's). A 
representative total flux spectrum from 2003 (bottom panel) is shown for comparison. Note the much narrower profiles of the Fe II and 
Ha emission lines in the polarized flux spectrum versus the total flux spectrum. 
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Fig. 9. — As in Figure Isl for the spectral regions surrounding the Na I D and Ha lines. Ro man numeral labels in the Na I D panel 
indicate the absorption systems described by IBoksenberg et ahl 119771) and IBoroson et alj Il99^ . 
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Fig. 10. — Aperture polarimctry of Mrk 231 on 1996 December 7 at B and R. In both colors there is a clear trend for P to decrease as the 
size of the measurement aperture is increased, indicating the inclusion of more unpolarizcd light from the host galaxy out to <~10" from the 
unresolved AGN. Triangles represent the polarization with a 3"xl0", east-west aperture that approximates a typical spectropolarimetric 
extraction aperture used for Mrk 231. 
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Fig. 11. — The polarized flux spectrum of Mrk 231 for 5000-6000 A in the observed frame (thin curves) convolved with gaussians ranging 
from no line broadening (the observed spectrum; bottom) to FWHM = 3000 km s _1 (top spectrum) in increments of 500 km s _1 . The 
thick curve is the observed total flux spectrum with the starlight contribution of the host galaxy within the observing aperture subtracted. 
Broadening the observed polarized flux spectrum by 1500—2000 km s _1 yields an emission-line sharpness similar to that of the total light 
spectrum. 



